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Abstract 

We derive the spectrum of the Dirac operator for the linear a- 
model with quarks in the large Nc approximation using renormaliza- 
tion group flow equations. For small eigenvalues, the Banks-Casher 
relation and the vanishing linear term are recovered. We calculate the 
coefficient of the next to leading term and investigate the spectrum 
beyond the low energy regime. 
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1 Introduction 



The QCD Dirac operator contains all information about the quark dynamics 
in QCD. Therefore its spectrum is of high interest. In the literature one 
considers the Euclidean Dirac operator which reads D = ^ + i g^. We will 
use Euclidean space time throughout this work. Since the Dirac operator 
is antihermitean, its eigenvalues are purely imaginary. With the eigenvalue 
equation 

D\^k) =i\k\^k) (1) 
one defines the spectral density 

p(A) = (5^5(A-A,)), (2) 

k 

where the averaging is over the gluon background. The first discovery about 
the spectral density has been made by Banks and Casher who connect 
the order parameter of spontaneous chiral symmetry breaking, i.e. the chiral 
condensate, to the QCD spectrum at zero eigenvalues. Since then, different 
approaches have led to a steady progress in this field. As the dynamics at low 
energies is mainly driven by chiral symmetry the behaviour of the spectral 
density in the vicinity of A = can be analysed using chiral perturbation the- 
ory. The first correction to the Banks-Casher relation was found by Smilga 
and Stern [0] 

o(X) - ^^^^"^^ + V ~ ^) lAI + O(X') (3) 

where V4 denotes the Euclidean 4-volume. The constant term is the Banks- 
Casher result, whereas the linear term vanishes in the case Nf = 2. 
Chiral random matrix theory considers the so called microscopic spectral 
density, in which all eigenvalues are rescaled according to the average spac- 
ing of small eigenvalues. It successfully describes the Dirac spectrum in the 
extreme infrared regime, where the spectral density is completely determined 
by the global symmetries of the Dirac operator (see e.g. [§]). Universal prop- 
erties of the microscopic spectral density have been derived and identified 
in lattice calculations 0, |^. Furthermore partially quenched chiral pertur- 
bation theory combines the results of chiral perturbation theory and chiral 
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random matrix theory 0. 

The Dirac spectrum for small eigenvalues is driven by chiral symmetry which 
determines hadron dynamics for small momenta. Thus chiral perturbation 
theory and chiral random matrix theory can be successful in describing the 
Dirac spectrum. 

The same holds for the Nambu-Jona-Lasinio (NJL) model which serves as 
an effective model for chiral symmetry. In this work we consider it in its 
bosonized form given by the linear cr-model. The QCD Dirac spectrum and 
the spectrum for the linear a-model which we analyze in this paper should 
coincide in the regime, where the linear cr-model is a good effective descrip- 
tion for QCD, that is for 

A < A . (4) 

Here A 1 GeV defines the momentum scale below which a hadronic de- 
scription of QCD may be useful. 

2 RG flow in the linear cr— model 

The Lagrangian of the linear cr-model in Euclidean spacetime is 

C = q[Z,^ + g{a + itt^^^)] q + {{d^af + {d^^f) + V{a, ^, m,) . (5) 

The fermion fields describe constituent quarks which appear in two flavours 
and three colours. The bosonic potential 

r(cT,7f,mg) = f/(<l>2) -ccr (6) 

contains a chiral symmetric part f/(<l'^) depending on the 0(4)-symmetric 
representation $ = (cr, vf) and a linear symmetry breaking term ccr. For 
later purposes we explicitly stress the dependence of the meson potential 
V{a, TT, TUq) on the current quark mass m^. At the ultraviolet scale, bosoniza- 
tion of the NJL model fixes the form of the meson potential VA(cr, yf, mq) to 
be 

V^{a,n,mq) = '^^'-ca. (7) 
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The linear symmetry breaking term c is related to the current quark mass 
ruq by 



m 



A 



m 



(8) 



We use renormalization group (RG) flow equations for the linear cr-model 
with quarks in the Schwinger proper time formalism. These flow equations 
have been discussed in various papers. For a general review on RG flow 
equations see and for the Schwinger proper time approach used here see 
[^]. The evolution equations provide a framework to study nonperturbative 
theories by integrating out smaller and smaller momentum shells above a 
limiting infrared cutoff k. This cutoff enters in the form of a smooth regulator 
function /fc(r) in the Schwinger proper time representation of the logarithm 
containing one-loop fluctuations. We consider the set of regulator functions 
given by 

f(-\re) = j2^-^e--'^\ (9) 

The flow equations which we present in this paper are obtained with n = 2. In 
order to check the cutoff dependence of our results we use regulator functions 
with different n. Generalization of the equations is straight forward. 
Note, the parameter c is scale independent, since a linear symmetry breaking 
term does not evolve under RG flow |^ . Therefore only the chiral symmetric 
part f/fc(<l>^) of the meson potential Vk{cr, vr, nig) = f/fc($^) — ca evolves under 
flow of the IR cutoff scale k. 

In this paper we restrict ourselves to the flow equations in the large Nc 
approximation, which have been considered recently [|TU| and whose results 



have proven equivalent to the standard selfconsistent NJL large approach 
pi] , [T^ . They are particularly simple and read 

dk 87r2 A;2 + $2 ' ^ ^ 

dZ^k _ NfN, k^ 
dk 47r2 (P + $2)3 • ^ 

The couplings Zq and g in eq. (|^) do not evolve in large Nc and have been 
set to 1 for convenience. 
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An additional flow equation which determines the chiral condensate can be 
obtained by differentiation of the partition function Z with respect to an 
appropriately introduced source term IS.qq which probes the chiral condensate 



(12) 

A=0 



This flow equation is derived analogously to the other flow equations by 
introduction of the Schwinger proper time cutoff in Z and reads 

dk Att^ + $2)2 • 

In large Nc all flow equations can be solved analytically. The initial conditions 
and solutions are discussed in |jlO[ . Denoting $ = {a, vf ) one finds for the 
potential (at k = 0) 

V{a, vf, m,) = V^{a, vr, m,) + ^ (^$^ log (^^^^) - $^A^) • (M) 

The minimum of this potential yields the vacuum expectation value $o of 
the bosonic field which is nonvanishing due to the spontaneous and explicit 
breaking of chiral symmetry. The physical value for the quark condensate is 

obtained as (gg)($o)- Its flow is shown in Fig. 1. 

1 

We obtain an infrared value {qq)fiG ~ 253.6 MeV, which compares well to the 

1 

value obtained from the Gell-Mann-Oakes-Renner relation {qq)GOR ~ 255.2 



MeV. This result is not surprising (cf. ||T2|) since the initial conditions of the 
flow equations have been adjusted to reproduce = 93 MeV and Mq = 320 
MeV at = using a current quark mass rriq = 7 MeV at the ultraviolet 
scale. 



3 Extracting the Dirac operator spectrum from 
the flow equations 

While the QCD-Dirac operator D is antihermitean, the Dirac operator in 
the linear a-model 

D = ^ + g{a + ifTr-f5) (15) 
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Figure 1: The evolution of the chiral quark condensate as a function of the 
renormahzation scale k with an explicit symmetry breaking due to a current 
quark mass of nig = 7 MeV. 




can have complex eigenvalues. In the flow equations we have neglected the 
imaginary part of the fermion determinant, which is related to anomalous 
processes [0. In this approximation the eigenvalues of the Dirac operator 
in the linear a-model are purely imaginary in coincidence with the QCD 
spectrum. For general complex eigenvalues z = x+iy it is useful to consider 
the two dimensional spectral density [0, [1^ 



(16) 



which can be obtained from the complex resolvent |[T^ 



G{z) 



Tr- 



z-D 



Tr— log [z — D 



d 

— log T)ei[z-D 



d 

— log / DqDqexp 



d'^x [q [ z — D) q 



(17) 
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The complex eigenvalue z is the generator in this resolvent, which in our low 
energy effective theory involves an average over the meson fields denoted by 
(...)$. In large approximation there are no meson fluctuations. Therefore 
the averaging over the meson fields can be interchanged with the logarithm 
in the last line and one finds 

G{z) = -V,l,qq){z) . (18) 

The connection between G(^z^ and P{x,y) can be derived with eq. (|l^) 

G{z) = 

X — x' + i{y — y') 

dx' dy' ^ . (19) 

X — x' + i[y — y') 

Combining eqs. (|TBp and (|T^) yields 

G(.. .) = -vm = . (20) 





This relation can be inverted |]T3| and gives 



where 



P{x,y) = ^ (21) 

vr oz 



dz* 2 \dx dy 

The complex eigenvalue z enters in the resolvent eq. (|r^) just as an explicit 
quark mass. Therefore no new flow equations have to be evaluated. After 
bosonization z appears in a linear symmetry breaking term in the meson 
potential via eq. (H) with niq = z. 

As the resulting meson potential Vfc(cr, vf, z) becomes complex, we must com- 
ment on the meaning of <l>o,fc(^)- Instead of restricting ourselves to positive 
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real values $o,A:, we now have to allow also complex fields ^o^k{z)- Instead of 
finding a real minimum of the potential Vk{(T, vr, rriq) we have to determine a 
complex saddle point, which minimizes the real part of the complex meson 
potential Vk{a, vr, z). 

From eq. ( ^1]) it is obvious, that the two dimensional spectral density is zero 
if {qq){z) is analytic in z. In presence of the complex quark mass z the chiral 
condensate at A; = becomes in the large limit with eq. 



{qq){z) 



N, 



c 



47]-2 



+ 2<i>o(^riog 1 + 



A2 + $o(z)2 uv , ^ ^ , ^^^^ 



(23) 



This expression is analytic if $o(-2) is analytic and Re {^q{z)) ^ 0. As $o(^) 
is the bare pion decay constant in presence of z and therefore minimizes 
Re (V^(o", vf, z)), the second condition Re($o(^)) 7^ is always true. But 
the global minimum of Re (V^(cr, vr, z)) can switch if z changes infinitesimally. 
This is exactly what happens and leads to a discontinuity in $o(-2^)- Fig 
demonstrates this behaviour. 

For z = it follows from eq. ([l4|) in the chiral limit that with $0 as a 
global minimum of Re (l^(cr, vr, z)) also — $0; '^'o ^"^^ ""^o global minima 
of Re vr, z)). Additionally they are saddle points of l^(cr, vf, 2;). For 
either z = x or z = iy the global minimum is still twice degenerate. In case 
of an infinitesimal x the sign of x chooses the minimum, i.e. Re ($0) changes 
its sign under the transformation from x = +e to x = —e. The same holds for 
infinitesimal y and for Im ($0)) but in this case an infinitesimal y only gives 
rise to an infinitesimal Im($), consequently no discontinuity arises in ^{z). 
We are thus left with the result, that P{x, y) can only be nonzero for x = 0. 
The position of our Dirac spectrum thus coincides with the QCD spectrum. 
Using the ansatz P{x,y) = 6{x)p{y) we find from eqs. (pOD, ( pT] ) and (|22| ) 



V4 _ 

piy) = TT- lim iw - e) - {qq) iw + e)] • (24) 

^vr e^O 

Now let us consider the flow equations for the potential and the chiral con- 
densate. In the large Nc limit they reveal a simple connection, namelyQ 

dk{qq)kmz)) = dkdMHz)) . (25) 

^ To be correct we should have written dk{qq)k{cr{z),Tf)\^^Q = 9^9^ ?7fe(cr(z), 7r)| -^p. As 
the expectation value of the pion field vanishes the notation above is unambiguous. 
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Figure 2: Projection of the real part of the meson potential depending on a. 
The sign of the explicit symmetry breaking term induces the choice of the 
global minimum and hence the vacuum expectation value. The solid curve is 
the potential without explicit symmetry breaking. The dashed and dashed 
dotted curves result from x = +e and x = —e respectively. 

2 

With the initial conditions {qq)A = and [/a = at = A this equation 

can be integrated out and yields 

{qq) ($(z)) = d^U mz)) - ml<!>{z) . (26) 

Eq. (|2^) gives the spectral density in terms of the discontinuity of the gener- 
alized chiral condensate {qq){z) across the real axis, which can be calculated 
from the saddle point (^o{z) of V{^,z) fulfilling d^V{^, z)\^^^^^ = 0. With 
eq. (^) we find in the chiral limit 

Piy) = ^ 1™ [*(^ y + e) - y - e)] 

= ^Rem^y + e)) . (27) 
vr 

The result is plotted in fig. (|l]). 

Next we derive the behaviour of the spectral density at small eigenvalues. 
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Figure 3: The spectral density p(A) of the Dirac operator D computed in the 
hnear cr-model with quarks. Note the suppressed zero on the scale for V^^p. 



Therefore we have to evaluate $o(^2/) for small y. The value ^o{iy) := $o(0) + 
a{y) + iP{y) is determined by the saddle point equation 



d^V{<^,y)\ 







— irrijy y = . 



(28) 



a{y) and [3{y) can be expanded in a Taylor series. The coefficients are ob- 
tained from 



dy^"-'^ 
We find 







j/=0 



^a{y) = $o(o) + i 



m 



A 



4'=*o(0) 



imidium 



n = 0,1,2. (29) 



*=*o(V + ... (30) 



$=$(0) 
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Figure 4: Behaviour of the spectral density at small eigenvalues, which ex- 
hibits a quadratic rise due to the absence of a linear contribution in the case 
Nj = 2. 



With eq. ( pT]) the spectral density 



p{y) 



IT 



Re[*„(0)l + i'"»"*^'*>l*-«™..= 



Slum 



1 3 



$=<I>o{0) 



+ 0{y') 



(31) 



can be calculated explicitly for small eigenvalues, see eq. (JT^) and ( p6D 

(32) 



UTT^ml (A2 + $0(0)2)*^ - Ar,A6 (SA^ - %{0f) (A^ + %{Off 



[An^ml (A2 + $o(0)2)' - A^.Ae] 

The constant coefficient is the Banks-Casher result |l|. Since our computa- 
tion is done for Nj = 2 there is no linear term in agreement with the findings 
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of Smilga and Stern for Nf = 2 0|. Numerically, the quadratic coefficient is 
87 MeV for the choice of parameters given in The quadratic term in y 
can be related to the correlation function of three scalar cr-quanta. In the 
numerator of eq. (^) one sees the effective coupling constant contained 
in our potential multiplied by the cube of the cr-propagator ^ at zero 
momenta. This result can also be found in the paper of Smilga and Stern 



4 The Dirac operator in the strong coupling 
regime of QCD 

We have presented a method to obtain the spectral density of the Dirac op- 
erator for the linear a-model in the large Nc limit. Because of asymptotic 
freedom in QCD, we expect p(A) oc for asymptotically large eigenvalues, 
since for a free theory the spectral density is only determined by the phase 
space. This situation of a non interacting theory can be tested within our 
approach in the limit m\ ^ oo. In this case the NJL four fermion interaction 
G oc l/m\ [|10| vanishes. As shown in fig. the resulting spectral density 
becomes the density of a free theory. For A > A the phase space is cut off 
sharply. 

If the interaction strength G is turned on by a finite mA one sees that quark- 
quark interactions lead to a level repulsion increasing the density for small A 
and pushing some strength above 1 GeV. Varying n and thereby the shape of 
the cutoff functions /^"^ in eq. (|^) we adjust the cutoff parameters m\ and 
A in such a way that f-jr and Mg are kept constant. This way we can test the 
sensitivity of the spectral density to the cutoff function. We see that below 
1 GeV the choice of the cutoff function leads to stable results within 10%. 
In the realistic case of QCD we expect that gluonic interactions above A 
influence the spectrum nearby and below. Therefore, the decrease of the 
spectral density p(A) above A ~ 0.8 GeV is probably an artifact of the model 
space, which is resticted by the soft cutoff. 

The interactions in the strong coupling regime of QCD induce level repul- 
sion. Higher modes are strongly suppressed, while the increase of the average 
eigenvalue density for smaller eigenvalues drives the chiral condensate which 
appears on the low energy end of the spectrum. It seems that level repulsion 
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Figure 5: The spectrum for different values of m\. As m\ approaches infinity 
the curve resembles more and more the free spectrum with a sharp cutoff at 
A. The spectra are plotted for mA = .228, .5, 1.5, 4, 10, 50, 2000 GeV. 
The free spectrum differs from the spectrum with mLambda = 2000 only for 
values of A very close to A. 

and condensation are the dominating effects for the Dirac operator spectral 
density in the strong coupling regime of QCD. 

5 Conclusions 

In this paper we have presented a new method to obtain the Dirac spectrum 
with renormalization group flow equations. We have used proper time flow 
equations for the linear a-model with quarks of two flavors. In the large 
Nc approximation we have found a simple connection between the meson 
potential and the chiral condensate which allowed us to calculate the 

Dirac spectrum directly from the full effective meson potential V{^, z). We 
have recovered the Banks-Casher relation for small eigenvalues. The next to 
leading order coefficient is 87 MeV. 
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Figure 6: The spectral density for different regulator functions fl^\ The 
uppermost curve corresponds to n=2. The others are respectively n = 
3, 5, 10, 15, where the last two are already nearly equal. 



The calculation of p(A) should be vahd up to several hundred MeV. Thus 
our approach allows to enter a new regime in the description of the Dirac 
spectrum which may complement other results from chiral perturbation the- 
ory and chiral random matrix theory. Compared to a free Dirac theory we 
have found the eigenvalue density suppressed for large eigenvalues and in- 
creased for small eigenvalues. All results are obtained in the infinite volume 
limit. Since the Schwinger proper time cutoff functions are simple to handle, 
the same calculation can be done for finite Euclidean space time volumes, 
where the evolution of the flow equations enters decisively. Such a calculation 
could then be compared to QCD lattice computations [|T5|. Further we plan 
to extend the calculation beyond the large Nc approximation. The renor- 
malization group flow equations of the linear a-model with quarks provide a 
framework to include meson loops in the NJL-model without an additional 
new cutoff [l^ . 
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